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Abstract 

Intriguing routes to chaos have been experimentally observed in semicon- 
ductor superlattices driven by an ac field. In this work, a theoretical model 
of time dependent transport in ac driven superlattices is numerically solved. 
In agreement with experiments, distorted Poincare maps in the quasiperiodic 
regime are found. They indicate the appearance of very complex attractors 
and routes to chaos as the amplitude of the AC signal increases. Distorted 
maps are caused by the discrete well-to- well jump motion of a domain wall 
during spiky high-frequency self-sustained oscillations of the current. 
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Nonlinear dynamics of weakly coupled semiconductor superlattices (SLs) driven by dc 
and ac biases has been the research topic of many experimental and theoretical works . 
The nonlinearity manifests itself in many physical situations as for instance, in their trans- 
port properties under finite DC or AC bias. Under appropriate DC voltage bias, the current 
through the SL displays natural oscillations caused by creation, motion and annihilation 
of domain walls (DWs) in the SL 0,^]. The observed oscillation frequencies range from 
kHz to GHz. Superimposed on a smooth current oscillation, there appear faster current 
spikes whose frequencies are typically in the GHz regime P]. Each spike reflects the well-to- 
well motion of DWs causing the self-oscillation, and therefore their frequency is related to 
the characteristic tunneling time. Theory and numerical simulations show that DW may be 
charge monopoles or dipoles although experimental evidence shows that DW are charge 
monopoles for self-oscillations observed in the available samples [^. AC driven SLs display 
a rich dynamical behavior including quasiperiodic and chaotic oscillations with nontrivial 
spatial structure pHTSl . 



Many studies fix the frequency of the AC drive as the golden mean number (1 + \/h)/2 ^ 
1.618 times the frequency of the natural oscillations (i.e., the frequency ratio is an irrational 
number hard to approximate by rational numbers), which is convenient to obtain complex 



dynamical behavior In this case, the system presents a rich power spectrum, a com- 
plex bifurcation diagram and different routes to chaos including quasiperiodicity, frequency 
locking, etc. |9Hr3|. First return or Poincare maps are used to analyze unambiguously the 
underlying attractors |l^,|T^. In the quasiperiodic case, Poincare maps usually consist of 
smooth loops, whereas they are a set of discrete points in the case of frequency locking. 
More exotic Poincare maps resembling distorted double loops in the quasiperiodic case have 
been experimentally observed in middle of the second plateau of the I-V characteristic of a 
SL |T2| , p!3| . At the onset of this plateau, Poincare maps are smooth and not distorted. The 
origin of distorted maps was not understood at the time of their observation, although dis- 



order and sample imperfections were invoked |T2[ . Luo et al [O showed that a combination 



of signals with different frequency was needed in order to reproduce experimentally observed 



distorted double layer Poincare maps. The origin of this combination was not ascertained 
in that work. 

The aim of the present work is to establish theoretically that high-frequency current 
spikes of the self-oscillations give rise to these exotic Poincare maps. In turn, current spikes 
are due to the well-to-well motion of DW during each period of the self-oscillations. Thus 
distorted Poincare maps reflect DW motion in AC driven SLs. 

Theoretical studies of nonlinear effects in weakly coupled SL are based in microscopic 
modelization of the sequential tunneling current between adjacent wells [^,0. Early models 
were discrete drift models with a phenomenological fitting of tunneling current and boundary 
conditions [Q. More sophisticated models have included the electrostatics at the contacts 
in a self-consistent framework . Depending on SL configuration and contact doping, we 
know that undriven self-oscillations are due to recycling of either monopole DWs or dipole 
waves [|8i. In the latter dipole wave consists of two DWs, one corresponding to 

a charge accumulation layer and another one corresponding to a depletion layer. Between 
these DWs there is a high-field region encompassing several wells, and the whole dipole wave 
travels through the SL. Recycling of such dipole waves gives rise to current self-oscillations 
similar to the well-known Gunn effect in bulk semiconductors [^. 

In this paper we have analyzed the tunneling current through a DC+AC biased SL by 
means of the microscopic model described in We have considered a 50-period SL 

consisting of 13.3 nm GaAs wells and 2.7 nm AlAs barriers, as described in Q. Doping in 
the wells and in the contacts are A^^ = 2 x 10^'^ cm^^ and = 2 x 10^^ cm^'^ respectively. 
With these doping values, self-oscillations are due to recycling of monopole DWs. We choose 
not to analyze the original 9/4 sample where the experiments were performed |T^. In this 



sample, the second plateau of the I-V characteristics ends at a resonant peak due to F — X 
tunneling. In the simpler SL configuration chosen here, the X point does not contribute 
to the peaks at the first and second plateaus. We analyze the sequential tunneling current 
1^ with an applied voltage, V(t) = VAcit) + Vac, where Vac(^) = Vac ^^^{^TcfAct)- The 
applied AC frequency /ac has been set to the golden mean times the natural frequency and 



is very small compared with typical energy scales of the system. Thus the AC potential 
modifies adiabatically the potential profile of the SL and our sequential tunneling model 
holds [gl. 

Fig. |1] shows the evolution of the current through the SL, its Fourier spectrum and its 
Poincare map for Voc = 4.2 V and Vac = 19 mV. These values correspond to the onset of 
the second plateau of the I-V characteristic curve. The current trace of Fig. |](a) is quasiperi- 
odic and does not present observable superimposed high-frequency oscillations (spikes). The 
natural oscillation near the onset of the plateau is caused by monopole recycling very close 
to the collector contact. Thus the DW does not move over many wells and the current trace 
does not present an appreciable number of spikes. In the power spectrum of Fig. |l](b) there 
are contributions coming from the fundamental frequency /o=39 MHz, the frequency of the 
applied AC field Jac, the combination of both and their higher harmonics. The Poincare 
map depicted in Fig. |I|(c) is a smooth loop with a nontrivial double layer structure indicat- 
ing quasiperiodic oscillations. If we increase the DC bias up to the middle of the second 
plateau, Vdc = 5.1 V, the undriven oscillation is due to recycling of monopole DWs which 
propagate periodically through part (approximately 40%) of the structure and disappear 
at the collector. A similar motion is observed in the AC driven case. The corresponding 
current trace may show high frequency spikes depending on the chosen initial field profile. 
A flat initial field profile gives a spiky current trace until the latter settles to the stable 
oscillation (over which there are no appreciable spikes); see the earlier part of Fig. H(a). 
The power spectrum and Poincare maps corresponding to this case and to the previous one 
(Figs. |l](b) and (c)) are markedly different. Fig. 0(b) shows that the power spectrum of the 
case with spikes has a greater harmonic content than that of the case without spikes (Fig. p, 
presenting many peaks corresponding to higher harmonics of the fundamental frequencies 
(i.e., the low fundamental frequency oscillations, /o = 33 MHz, and the spikes, 500 MHz), 
the applied frequency and combinations thereof. The corresponding Poincare map consists 
of a distorted loop with a double layer structure which shows a strong similarity with pre- 
vious experiments [|T2i|l3| ; see Fig. ^(c). From these numerical observations, we conclude 



that distorted Poincare maps are linked to spiky current traces, even if such traces change 
to smooth ones after a transient. 

The previous conclusion may be reinforced if we change the emitter doping so that 
undriven self-oscillations are caused by dipole waves, for which the current traces are more 
spiky 1^. Thus we carried out numerical simulations with a smaller contact doping, N^. = 
2 X 10^^ cm~^. The frequency of the natural oscillation is now reduced to 4 MHz and 
the AC intensity is 2 mV. The current at the middle of the first plateau (Vdc=1-5 V) is 
shown in Fig. ^(a). At finite time, I{t) presents dipole-like oscillations and superimposed 
finite amplitude spikes. The Poincare map. Fig. 0(c), is much more complicated than in the 
previous case, showing three well defined distorted loops. Since loops in the Poincare map 
are due to combination of strong enough signals of different frequencies [0, the greater 



strength of the high-frequency spikes gives rise to the additional loop structure and higher 
harmonic content (Fig. |^(b)). 

As we mentioned above, the double layer structure of Poincare maps indicates nontrivial 
attractors for the quasiperiodic case. We have calculated the multifractal dimensions of 
the attractors, Dg, for the three cases discussed above; see Fig. ^. In all cases, they corre- 
spond to strange attractors [|TU],|TR|I7| whose Dg presents the knee-like structure typical of 



chaotic attractors with multifractal dimensions |T^. A detailed discussion of the attractor 
dimensions will be presented elsewhere. 

In summary, we have analyzed theoretically the time dependent current through a AC 
driven SL. We have characterized intriguing Poincare maps of the quasiperiodic oscillations. 
The strange attractors which define these Poincare maps have their physical origin in the 
complex dynamics of the domain wall. We have shown that distorted loops in the Poincare 
maps are related to the presence of high frequency spikes in the current traces. Their fre- 
quencies combine with the AC frequency and the low natural frequency to yield a richer 
power spectrum. Since spikes are associated to extended motion of the DW, we can under- 
stand why distorted Poincare maps are observed at the middle of a I-V plateau and not 
at its beginning, where the monopole moves over a too small part of the SL. The case of 



natural self-oscillations due to dipole recycling is different. There the DW of the dipole are 
generated at the emitter contact and move over the whole SL. In such a case, spikes are 
more prominent than in the monopole case, and a higher harmonic content and distorted 
Poincare maps appear. Our results may help explaining recent experimental evidence show- 
ing complex Poincare maps and intriguing routes to chaos 
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FIGURES 
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FIG. 1. (a) I(t) for Vdc = 4.2 V, /o = 39 MHz, Vac = 19 mV. Spikes are not resolved, (b) 
Power spectrum. Notice that higher harmonics of the fundamental frequency are barely formed, 
(c) Poincare map, constructed by plotting the current at the (n+l)st AC period versus the current 
in the preceding period. 
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FIG. 2. (a) I(t) for Vdc = 5.1 V, /o = 33 MHz. Spikes are present in the transient regime 
(see inset), (b) Power spectrum. Some higher harmonics of /o can be observed, (c) Poincare map. 
Notice that the loop is somewhat distorted due to the presence of higher harmonics in (b) . 
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FIG. 3. (a) I(t) for Vdc = 1-5 V, /o = 4 MHz, Vac = 2 mV. Spikes are superimposed on the 
current throughout the signal (see inset), (b) Power spectrum. Higher harmonics of /o contribute 
with a finite amphtude to the power spectrum, (c) Poincare map. The distortion is greater than 
in Fig. 0(c) (see text). 
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FIG. 4. Multifractal dimension. The labels (1), (2) and (3) correspond to the attractors in 
Figs. 0(c), Q(c) and ^(c), respectively. 
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